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We demonstrate the surprising integrability of the classical Hamiltonian associated to a spin 1/2 
system under periodic external fields. The one-qubit rotations generated by the dynamical evolution 
is, on the one hand, close to that of the rotating wave approximation (RWA), on the other hand 
to two different "average" systems, according to whether a certain parameter is small or large. 
Of particular independent interest is the fact that both the RWA and the averaging theorem are 
seen to hold well beyond their expected region of validity. Finally we determine conditions for the 
realization of the quantum NOT operation by means of classical stroboscopic maps. 
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I. INTRODUCTION 



Advanced NMR techniques arc able to manipulate 
qubits in order to implement several quantum logic oper- 
ations. In particular, by means of radio- frequency mag- 
netic fields, nuclear spin 1/2 can be controlled to the 
realization of the quantum NOT-gate, one of the funda- 
mental logic block in quantum computation On the 
other hand, in the field of quantum optics, the advent of 
strong laser pulses brings up the possibility of investigat- 
ing the fundamentals of the radiation-matter interaction. 
In this scenario, models like a two-level quantum system 
interacting with an external classical field have achieved 
a paradigmatic place, since it has been succefully applied 
by the theorists to make predictions in NMR quan- 
tum computation and quantum optics 0. 

From a mathematical point of view its importance 
emerges mainly from the low dimensionality of its Hilbcrt 
space, which allows analytical analysis of the dynamics in 
many cases of interest. For periodic external fields, e.g., 
such systems have been solved analytically only in the 
rotating wave approximation (RWA) , which restricts the 
problem to both resonance and weak coupling regimes. 
Far off this situation, sophisticated perturbative methods 
have been applied 0, 0, |(| and exact solutions have been 
found only for specific periodic fields 0. An alternative 
approach for studying exactly the two-level problem is 
provided by the picture of the classical gyromagnet 
The dynamics of the matrix elements of a pure density 
operator is mapped on a classical Hamiltonian and meth- 
ods of dynamical systems may be employed to study the 
dynamics of the system. Alternatively, it is possible to 
define a classical dynamics over the Bloch sphere @ . 

In this work we focus on these classical frameworks to 



study relevant questions as the integrability of the under- 
lying classical dynamics, the validity of the RWA (as well 
as its derivation from the averaging method) and regimes 
for the realization of a unitary NOT operation. The pa- 
per is organized as follows. In section ^ we define the 
two classical frameworks on which the analysis will take 
place. In section Hill we demonstrate the integrability of 
the classical dynamics associated to unitary two-level sys- 
tems under arbitrary external fields. In sections IIVI and 
[V]we present several results concerning the phase space 
aspect of two important fields commonly used in exper- 
iments. In section IVII we use the results offered by the 
stroboscopic maps to make predictions about the quan- 
tum NOT operation. Section IV11I is reserved for some 
concluding remarks. 



II. THE CLASSICAL FRAMEWORK 

The dynamics of a spin 1/2 system interacting with an 
external time-dependent magnetic field B(f) in the dipole 
approximation is given by the Hamiltonian 



H(t) 



-B(t)-E, 



(1) 



in which X = (<7i, a%, 03) is the vector composed by the 
Pauli matrices. For simplicity we have adopted H = 1, 
so that magnetic field is given in units of frequency. The 
dynamics of the density operator p(t) = \tfj(t))(ip(t)\ is 
governed by the von Neumann equation, ip = [H(tYp\. 

As mentioned, according to Feynman et al @|, it 
is possible to formulate the dynamics in terms of a 
classical Hamiltonian systems. Consider the following 
parametrization for the quantum dynamics in terms of a 
classical vector (So, S): 
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iS 2 
-S 3 
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in which 1 stands for the unity matrix. The trace (Tip = 
1) and purity (p 2 = p) conditions impose that Sq = 1 
and S 2 = S'f + S 2 + S 2 = 1, respectively. The quantum 
equation of motion for p is then transformed in a system 
of differential equations described by 



7/7 



= S x B, 



(3) 



which may be obtained from the classical Hamiltonian 



«(*) = -B(t) • S. 



(4) 



Equations and J3J define the classical geometric pic- 
ture associated to quantum two-level systems: the unit 
vector S precesses around the vector B just like a classi- 
cal gyromagnet precesses in a magnetic field. 

Following reference we consider the unit sphere <S 2 
with the usual angular coordinates ^ 6* ^ 7r, ^ ip ^ 
2ir, and let S = (sin 9 cos tp, sin 9 sin (p, cos 9) . Introducing 
p = ip and q = — cos 9 as canonically conjugate variables, 
we may write 



S = (yl - q 2 cosp, yl - q 2 ship, — g), 



(5) 



with the usual Poisson brackets {Si, S 2 } = S3 (plus ciclic 
permutations). By (@J) and JSJ we finally write 



H = - [B x {t) cosp + B 2 (t) ship] ^l-q 2 - B 3 (t)q, (6) 



with equations of motion given by 



q = [Bi ship - B 2 cosp] yl - q 2 , 
p = — \B\ cosp + B 2 sin p] 



(7a) 
(7b) 



in which the amplitudes Bi are functions of time. 

In our analysis it will be useful to express the quan- 
tum states in terms of the classical canonical pair (q, p) . 
Consider a general quantum state of a single qubit as 
= a|+)+&e~'"^| — ) being a and b positive real numbers 
satisfying a 2 + b 2 = 1 and <f> a the relative phase. Com- 
posing the corresponding density operator, p = 
and comparing with J2J) and (JSJ, one may write 



1 — Q, 



1 



9 



(8) 



Thus, correspondence between classical phase space and 
quantum states becomes immediate. For instance, for 
q = =pl, we have S3 = ±1 and \ip) = |±), with obvi- 
ous interpretation. Notice that the relative phase of the 
quantum state has the status of a conjugate momentum 
in classical phase space. It indeed emphasizes its dynam- 
ical relevance. 

We will also consider a perpendicular state \ip±), de- 
fined by (ip±\^p) = 0. According to JSJl it reads 



1 



1 -g ,( 

— ^ — e 



-*>!->. 



(9) 



Thus, the transformation 

1^) \i>±), (io) 

implies in the following classical counterpart 

(p,q)^ {p + TT,-q). (11) 

By © we see that such a condition of orthogonality im- 
plies that S — > — S. These considerations will be essential 
for our analysis of the quantum NOT operation in two- 
level systems. 

III. INTEGRABILITY 

The classical Hamiltonian @ may be written in a two- 
dimensional autonomous form by means of Howland's 
method [Toj . In this case, since the new Hamiltonian is 
composed by an integrable part added to a non-integrablc 
perturbation, chaotic behavior is expected to be present 
Q. However, as we show now, the unitarity of the un- 
derlying quantum dynamics guarantees that this is not 
the case. Consider the following distance in the Bloch 
sphere: 



D(t) = \\Si(t)-S 2 (t)\\, 



(12) 



where different indices refer to different initial conditions. 
By © it is easy to show that D may also be written as 



D(t) = v /2Tr[ (0l (<)-p 2 (t)] 2 , 



(13) 



where different indices refer to different initial states 
Pi, 2(0). Since the dynamics is unitary, we may write 



pi,2(t) = U(t)p lt2 (0) U\t), 



(14) 



in which U is the propagator satisfying 1U = H(t)U(t). 
Equation (|12|) allows us to define a Lyapunov exponent 



as 



A 



lim lim — In 

D(0)^0 t^oo t 



D(t) 



D(0)\ ' 



(15) 



which measures the mean exponential departure be- 
tween two arbitrarily close initial conditions on the Bloch 
sphere. By (|14|) and (|13fl it is immediate that D(t) = 
D(0) and thus, by JTSJ, A = 0. When a full set (i.e, 
equal to the number of degrees of freedom) of constants 
of the motion in involution does not exist, so that we 
have a system which is not integrable, sensivity to ini- 
tial conditions is certain to exist at least in some region 
pf the phase space, i.e., there is inevitably chaotic be- 
havior. Thus we may regard the proof that A = as a 
proof of integrability, although it may be difficult to find 
explicitly an additional constant of the motion for the 
eneral autonomous system equivalent to (0} (see, e.g., 
, for the description of these autonomous systems). It 
is worth emphasizing that our proof of integrability is 
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valid for an arbitrary time-dcpcndcnt field, including the 
quasiperiodic case. This contradicts the basic assertion 
of and confirms the result of 0] . More importantly, 
we show the real reason for the integrability observed in 
using double Poincare sections. Notice also that a 
field in Q depending of N incommensurate frequencies, 
would yield, by Howland's method ^(Jj a A^+l degrees of 
freedom autonomous system. Our demonstration asserts 
that such a system is integrablc for all N, a surprising 
result. 

We close this section with some remarks which may 
help to connect the present model with the fundamental 
model of interaction between a (two-level) atom and the 
quantized electromagnetic field in the dipole approxima- 
tion |l3j , which we shall take as one- mode for simplicity: 

H = ua'a <8> 1 + -y 1 <g> a z +g (a + a f ) ® a x , (16) 

where a and a 1 are the annihilation and creation opera- 
tors, respectively, and 1 is the identity matrix (properly 
defined in each subspace). Although "quantum integra- 
bility" is not a well-defined concept in general, two natu- 
ral definitions, with direct analogue in classical mechanics 
are possible: 1) l|16|) is integrable if there exists a uni- 
tary transformation U\ such that U^HU^ = H(N,a z ), 
where N = a^a and a z are two "quantum actions" (com- 
muting operators, whose spectrum and eigenvectors are 
known explicitly); 2) l|16f) is integrable if there exists a 
unitary transformation U2 such that the spin and bosonic 
degrees of freedom decouple. It has been proved in Q 
that (|Tfi|l is integrable according to 1) for both large and 
small coupling g. In reference |l5j it has been remarked 
that 



H = 



A- 



U 2 ] HU 2 



= — a z cos (ko) a) + oj a^a + g (aJ + a) , (17) 



where 



U-2 = exp {— itt (a x — 1) a^a/2}. 



(18) 



Thus, by (|17f) . in each sector corresponding to the eigen- 
values ±1 of a z the system is equivalent to a one- 
dimensional model (equivalent to one degree of freedom 
in classical mechanics), and thus the spin and bosonic 
degrees of freedom are decoupled. Notice, however, that 
the resulting model is highly "nonsoluble" ! This fact re- 
flects upon the semiclassical limit of <|16[l . in which the 
"mean photon number" (properly defined) tends to in- 
finity, keeping the photon density fixed [l^ : it yields our 
model Q , with B = B mr given by (fT^TTf) , whose dynamics 
is known, since the seminal papers of Bloch and Siegert 
|l7j and Autler and Townes |l8j ] , to be highly nontrivial. 
For the system l|16|) there exist two commuting operators: 
-ff, given by (|16f> . and the parity operator 



II = exp < vk 



1 



a T a + - {a z + 1) 



The fact that II corresponds to a discrete symmetry 
without classical analogue might have suggested that the 
semiclassical limit Q, with (jSOJ, is not integrable, which 
we have seen not to be the case. 

Above, we have considered one atom only. A different 
limit which may be performed on (|16J) is the "many-atom 
limit" (properly defined [l9j|). The method of justi- 
fies rigorously the results of (2^, I2H |2^ . whereby the 
resulting coupled classical equations may be shown to 
display chaotic behaviour. 



IV. THE ROTATING FIELD 

In this section we will focus our analysis in a special 
case of periodic external field: the radio-frequency field 
0. This field, which we will call rotating (R) in virtue 
of a rotation symmetry around the direction 3, is defined 

as 

B R (t) = -2 (b q cos (u>t + (j>), B Q sin (art + B 3 ), (20) 

being <j> the phase of the field and Bq and B3 constant 
amplitudes. By © and l|20(l . we write the corresponding 
Hamiltonian as 



H(q,p,t) = 2B Q ^l-q 2 cos (p + - ut) - 2B 3 q. (21) 

The integrability of this model may be explicitly verified 
by means of a canonical transformation that produces an 
autonomous one-dimensional Hamiltonian. Consider the 
following generating function and respective transforma- 
tion equations pjj : 



F 3 (p,Q,t) = -(jp + 4>)Q + e{t)Q, (22a) 
dF 3 
dp 

dQ 
dF 3 



<1 - 

P = 

k = n 



Q, (22b) 
■p+(t>-6, (22c) 
H + BQ, (22d) 



being 8(t) = cut in this case. It is just the classical coun- 
terpart of the usual unitary quantum rotation given by 
exp (— iujto-3). Equations above yield 



K{Q,P) = 2B a ^/l - Q 2 cosP- 2f2Q, 



(23) 



(19) 



with Q = B 3 — w. Now we have an integrable Hamilto- 
nian, since /C is a constant of motion in a one-dimensional 
system. Therefore, the system R defined by l|21|) is man- 
ifestly integrable. Nonetheless, because the phase space 
(p, q) contains separatrices (as we will see latter) finding 
analytical solutions for the equations of motion may not 
be a simple task. The simplest (and more reliable) way of 
getting the solutions is solving the problem in the original 
quantum system and then implementing the correspond- 
ing transformations. We present the analytical solutions 
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in the appendix [X] In the next sections we will deal 
with the explicit solutions. For now it is sufficient for 
our purposes to understand the general characteristics of 
this integrable system. 

Before constructing the stroboscopic maps, we point 
out one further interesting information. By (|22ll : 



or 



H = IC{q ,po) - uq, 

dH _ 
dq 



(24) 



(25) 



which indicates that the time-dependent energy. 7i, keeps 
a linear relation with q for all times. 



A. Stroboscopic Map 

In the general case of a one-dimensional Hamiltonian 
system driven by an external periodic field, for the which 
analytical analysis is not possible, stroboscopic maps are 
useful tools to attest the integrability. In our case, al- 
though we have already proved the integrability, it will 
be useful constructing the stroboscopic maps in order to 
obtain some sights about the whole aspect of the inte- 
grable tori. As we will see, it will be essential in our 
analysis of the NOT operation. 

The stroboscopic map is obtained by marking the po- 
sition and its conjugate momentum at each multiple of 
the Hamiltonian period, T, which is defined by the re- 
lation H.(q,p,t + T) = 7i(q,p,t). A rapid inspection of 
the Hamiltonian J2TJ) yields T = ^. We then define the 
stroboscopic instants as = kT, for k 6 N. 

Since we have obtained the analytical solution for the 
system, the construction of the stroboscopic map is triv- 
ial. However, another procedure based on qualitative de- 
scription of the stroboscopic map will be more appropri- 
ate for the comparative analysis that will be realized in 
the next sections. Re-writing <|22d(l at the stroboscopic 
instants one obtains 



K, = 2B 0V /l-ft 2 cos ( Pk + cp)-2[Bz- - J q k , (26) 

in which denotes r(t/ £ ). This formula allows one to 
construct the stroboscopic map by means of the contour 
curves. It is an alternative indication of the integrability 
of the system, since the map is composed only by tori 
(contour curves). This observation avoids (numerical) 
dynamical evolution, although it does not offer informa- 
tion about the comensurability pattern of the torus (ra- 
tional or irrational), which may be obtained only by the 
analytical solutions. In fact, by means of the analytical 
equations of motion (see appendix we found the con- 
ditions of comensurability, given by the ratio — , being 
B = 2^/Bq 1 + ft, 2 the amplitude of the magnetic field in 
the rotating frame (see appendix[BJ) . Notice that rational 
and irrational tori cannot coexist in the same map, since 



the comensurability parameter, ^, does not depend on 
the initial conditions. Figure ^ illustrates these informa- 
tions by some numerical examples. 
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FIG. 1: Stroboscopic maps obtained numerically from the 
whole evolution (a,c,e) and respective contour curves obtained 
by {HI (b,d,f) for the R field with: B = 1, B 3 = 0, ut = 1 
and § = ^5 (a,b), B = B 3 = w = 0.5 and § = VE (c,d) 
and Bo = 1, B3 = -| , w = 89 and S = JL (e,f). In (e) we see 
a case of rational tori. d> = in all cases. 



B. Periodic Orbits 

At stroboscopic instants both set of variables, (p,q) 
and (P, Q) yield the same map, since Qk = q^ and Pk 
(mod 2ir) — p^. (mod 2ir). Then we can calculate the 
fixed points (periodic orbits of period one), denoted by 
(P, Q), by imposing that Q = P = in the equations of 
Hamilton. This yields 



P = In 



(I = 0,1, 2), 



Q ■ 



2Q 



(27a) 



(27b) 
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According to (jSJ the corresponding quantum state (in 
the rotating frame) is 

M - ± ^±p|->. (28) 

Putting (P,Q) in the Hamiltonian (|23|1 . we obtain the 
contour value K, — E, to which the fixed point belongs: 

E = ±B. (29) 

Equations (|28|l and i|29|) given the cigenstates and the 
eigenvalues of the quantum Hamiltonian associated to 
K,. Actually the eigenvalues arc ±B/2, as one may guess 
by (JTJ and J3J . This is a very interesting example of the 
correspondence between eigenstates and periodic orbits. 
Another interesting point is that the frequency f2 = B% — 
tt controls the vertical localization of the fixed point and, 
consequently, the symmetry of the map. For instance, in 
the resonance (f2 = 0) one obtains Q = 0, as can be seen 
in Fig. [2 

Finally we point out that initial conditions like l|27|) are 
fixed point for all times in the phase space associated to 
variables (P,Q). However, this is not so for the original 
variables, since that by l)22ll we have (p(t),q(t)) = (P — 
(f) + uit, Q±). In fact, in these phase space the fixed point 
occurs only at stroboscopic instants. 



V. THE NONROTATING FIELD 

We are now interested in analyzing the nonrotating 
(NR) field defined by 

B NR =-2{B ,0,B 3 cosut), (30) 

being Bq and B% constant amplitudes. Actually this is 
one of the most common fields used in NMR experiments, 
mainly for quantum computation purposes , even more 
because fields like (|2(Jfl are difficult to manufacturate. 
The corresponding Hamiltonian reads 



H(q,p,t) = 2B0V 1 - I 2 cosp-2B 3 qcosu;t. (31) 

The absence of a rotation symmetry destroys the mani- 
fest integrability of the model. In the R system, Br = 
\\Br\ \ is constant, so that the time variation of the exter- 
nal field is given only by changes in its direction. Then 
one may choose a frame that rotates with the field fre- 
quency, getting a constant energy. Because in the NR 
system the field varies also in modulus it is not possible 
to eliminate time-dependence just by defining a rotating 
frame. In fact we are not able to find a canonical trans- 
formation that puts the Hamiltonian in an autonomous 
form. The quantum version of this system has been ana- 
lyzed in different contexts, but always in a perturbative 
way 



A. Numerical Analysis 

In this section we present some numerical results con- 
cerning the features of the stroboscopic map and its qual- 
itatively ressemblance with those of the R system. Con- 
sequently, integrability is numerically verified. 

The first interesting information that emerges from the 
numerical results is the linear relation between the time- 
dependent energy, H(q(t),p(t),t), and the coordinate q(t) 
at the stroboscopic instants tk = ^j-k. In Fig. [21 we illus- 
trate this information by means of stroboscopic bullets. 
It is also possible to notice, by the continuous flux, the 
complexity of the NR dynamics (compare to l24fl . 
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FIG. 2: Parametric plot between the time-dependent energy 
Tt(t) and the coordinate q(t) for the continuous flux (solid 
line) and for the stroboscopic instants tk (bullets). The pa- 
rameters used in this calculation were Bq = 1.0, B3 = 1.5 
and to — 3.0. 
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FIG. 3: "Stroboscopic parametric plot" of TL and q (sym- 
bols) for several initial conditions and respective fittings (solid 
lines). The slopes of the fitting curves (7 = 4.9559) are identi- 
cal and independent of the initial condition. The parameters 
used in this calculation were Bo = 1-0, B3 = 1.5 and w = 3.0. 

In Fig. |21we show the same type of stroboscopic para- 
metric plot for several initial conditions under the same 
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parameter interaction. It is possible to see that the slope 
(regression coefficient), 7, of each straight line has always 
the same value, i.e, it does not depend on the initial con- 
dition. On the other hand, the regression constant, £, 
does depend on the initial condition. 

These numerical results allow us to write = —7 
and 

Hk = £(qo,Po) - jqk, (32) 

in which AH = H k+ i - Hk, Aq = q k+l - q k and Hk = 
H(qk,Pk,tk)- Using 1)3 1|> we re-write equation above as 

S(q ,p ) = 2B 0V /f -qk 2 cos Pk - 2 (b 3 - |) q k , (33) 

which is the equation of the stroboscopic map for the 
NR system as a function of the numerical parameter 7. 
Comparison between <|32|) and (|24|l emphasizes the re- 
semblance between NR and R dynamics, at least at the 
stroboscopic instants. 

On the other hand, equation (|33|) indicates that the 
maps of the NR system must be "identical in shape" to 
those of the R case, as can be seen by comparing <|33[) 
and (|26(1 (with = 0). The difference is only in the reso- 
nance condition, which is now given by -B3 — ^ , instead of 
B3 — Tf. 7 is a frequency determined numerically by the 
fittings shown in Fig. |3J whereas u> is the own frequency 
of the applied field. The qualitative agreement between 
NR and R maps are shown in Fig. 0] which was made 
as follows: we determined 7 numerically by fittings in 
the NR system and then we used a R field with ui = 7. 
Another interesting difference can be noticed in these nu- 



(a) (b) 




FIG. 4: (a) Stroboscopic map for the NR system with Bo = 
1.0, Bi = 1.5 and u — 3.0 and (b) the corresponding map for 
the the R system with u> — 7 = 4.9559 (see eq. 1261 . This 
value was obtained by the fitting shown in Fig. [3] 

merical results: in the NR system, rational and irrational 
tori are allowed to coexist. 

The above results explain the fact observed in £j that 
the RWA accounts for several results qualitatively even 
outside the regime where it is supposed to be a good 



approximation [l3| . namely: i.) resonance and ii.) weak 
coupling. In particular the solution found in Q in the 
case of resonance has a frequency which is very close to 
the Rabi frequency. 

In the next sections we demonstrate the explicit in- 
tcgrability (i.e., the possibility of obtaining explicit ana- 
lytical solutions) of the NR system at some regimes of 
interest. 



B. Analytical Analysis: Average Dynamics 

The contour pattern (|33|l chosen by the NR system 
is not obvious. In fact it can not be inferred directly 
from the Hamiltonian (|31fl . The choice by one specific 
contour in the set of possibilities offered by the time- 
dependence of the system suggests the existence of an 
average dynamics. In this section we try to obtain some 
analytical evidence of such an idea. 

Consider the following arbitrary Hamiltonian 

H(q,p,t) = H (q,p) + f(u>t)V(q,p), (34) 

in which f{tut) is a periodic function with period T = — . 
We now calculate the time-derivative of the Hamiltonian, 

H = Vd t f = d t (Vf) - fd t V, (35) 

and then integrate the result over the arbitrary interval 
[^1^2] to get a time-dependent function given by 

H 2 -Hi = I2V2 - fiVx - / f{wt)V{t)dt, (36) 

in which we used the notation A{ = A{ti), with A stand- 
ing for H, V and /. Now, assuming that V(t) is not 
a constant of motion, we multiply and divide the last 
term of the equation above by V2 — V\ — J t 2 V(t)dt and 
reorganize the equation in the form 

n 2 -n x = (/ 2 - (f) 2 ,i)v 2 - (a - (/> 2 ,i)vi (37) 

in which we have defined the average 

J f(ut)V(t)dt 

</>a,i = H ■ ( 38 ) 

JV(t)dt 
ti 

This equation can be re-write as 

J f{ujt)V(t)dt = (/) 2a (v 2 - Fr) , (39) 
ti 

which presumes that the integral of fV may be identified 
to the product of the integral of V alone multiplied by a 
certain average. 
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At arbitrary stroboscopic instants, namely t\ = iT and 
< 2 = jT, with i and j G N, ||38J and J2SJl may be sufH- 
cicntly manipulated to yield 



</>i,< = 



and 



T 3-1 . 

J f(ujt) E V(t + nT) dt 

n=i 

Tj-l . 

/ E v(* + nT ) * 

ri=i 



(40) 



(Vj - V.) {f)j,i = (Vn+1 - V n ) (f) n+ l,„. (41) 



Convenient transformations of variable were performed 
in order to put the formulas in terms of integrals over 
only one period. 

Consider the behavior of the mean 14U|) at consecutive 
stroboscopic instants for long times (k 3> 1)- In such 
case we have 



k+l,k 



J f(ut)V(t + kT)dt 
o 

T 

J V(t + kT)dt 



(42) 



For any value of T there exist a value of k large enough 
to guarantee that kT ^> t. It allows us to expand V (kT+ 
t) in Taylor's series, being t the small parameter. This 
procedure leads to 



(fh 



+i,fc 



E 

n=0 



o 



E 



a?,v(t k ) 



1 (2 2L \ n + 1 
n+1 Ui 



(43) 



in which ip = cut. In this result all dependence with ui 
was made explicit. Then, the regime of high frequency 
becomes evident for well-defined functions / and V: 



Jim </) fe+1 , fe = i f f{ut)dt = f. 

(<=>!) J 



(44) 



In such situation, according to 1)41(1 . we have = /. 

Then, choosing t% = tk, ii = and putting results above 
in i|!r7|) we obtain 



n k = 



KoiqoiPo) + fV(q ,po) 



/o - / 



£(q ,Po) = n k - fo-fV k 



V k , (45) 



(46) 



with £{q Q ,pa) = 7io(9o,Po) + fV(q ,Po)- This result at- 
tests the intcgrability of the system, since it provides the 



analytical expression for the stroboscopic map. Notice 
that 



AH_ - 



(47) 



and thus (|45|) and l|47|l can be directly compared to re- 
spective results (|24|l and (|25|) for the integrable system 
R. Once the system is notably integrable for k ^> 1 it 
is expected that it is also true for every k. It must be 
so to guarantee that the map is composed by the same 
contours at any stroboscopic instant. 

The regime of low frequency can not be predicted di- 
rectly from the formulas developed in this section and nu- 
merical work is necessary. However, the high frequency 
regime illustrates the necessary and sufficient condition 
for the integrability of the system l|34|) : 



fe+i,fe 



= (/)• 



(48) 



This condition requires that the average defined by l|42|l 
be time-independent. One further observation should be 
done in order to emphasize that the average dynamics is 
weighted up by the time-derivative of the potential and 
not only by the time-average of the function /. This fact 
credits a non-trivial character to the identification of such 
an average. 

In the next section, we obtain some explicit regime of 
integrability by applying the average we defined above, 
the averaging theorem and the RWA. 



C. Analytical Solutions in the NR System 

1. High Frequency Regime 

In this section we consider the regime of high fre- 
quency, in which the frequency of the external field w 
is much larger than Bq or B3. In fact, we will work in a 
regime that allows us to disregard higher order terms in 

63- 

At first we investigate the terms 



1 



Jo 



2jv 



and 



B n = 



1 /2?r 



n + 1 V u) 



n+l 



(49a) 



(49b) 



that appear in the brackets in For f(ip) = cos<^ the 
integral above can be performed analytically, leading to 
a hyper-geometric function with a maximum at n — 2. 
In Fig. [SJwe show the behavior of w 3 A„/47r as a function 
of n for two values of u>. For high frequency regimes we 
indeed verify that only n = 2 contributes significantly to 
the sum, i.e., only Ai is non-null, though small. On the 
other hand, for B ni the major contribution comes from 
n = 0. as it is evident from its own definition. 
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FIG. 5: " 4 ^" as a function of integers n for uj = 10 (dashed 
line) and ui = 100 (solid line), in arbitrary units. The maxi- 
mum value of the function plotted does not depend on u>. 



According to Fig. |3J we can write A n = ^-S n .2- Then, 
(|43|l assumes the form 



(f] 1 dlV{t k ) 



(50) 



Now, by noticing that V = —2B 3 q and using the equa- 
tions of motion of the NR system, we get 



(fh 



k+l,k 



, = -4 



So 5 



(/>, 



(51) 



which is notably fc-independent. This result explicits the 
condition of high frequency that we supposed initially. 

For large values of uj (not infinity) we can replace / in 
(|45|l by the result Q51fl. Then, by comparing l|45[) with 
(I32fl . wc establish the connection between the analytical 
((/)) and numerical (7) results as follows 



7 = 2B 3 (!-</» 



(52) 



In Fig. |HJ the validity of this relation for high frequen- 
cies is verified, showing that our approximations are ad- 
equate. 

The analytical expression for the stroboscopic map 
in the regime of high frequency can be constructed by 
putting Vk = —2B 3 q k in (|45|l with the explicit form of 

7~tk'- 



£(q ,p ) = 2B oy /l-q k *cosp k - 2B 3 (f)q k , (53) 

with (/) given by (JHU). Notice that 7 -> 2B 3 and (/) -> 
in the limit of uj — > 00. It shows that the map will have 
periodic orbits localized in the line q = 0. In this limit, 
accordingly to both ffity and , the stroboscopic map 
reads 



£(<?o, Po) = 2B a/! - qk 2 cosp k - 



(54) 



The approach adopted in this section allows us to pre- 
dict the form of the stroboscopic map in the regime of 



FIG. 6: . Comparison between the numerical 7 and analytical 
2Bz(l — (/}) slopes for several values of uj. 



high frequency, defined precisely by H51|) . However this 
cannot be achieved for any periodic function /. In fact, 
it is possible to show that the average is fc-dependent for 
f(t) = sinujt. 



2. Strong Coupling Regime 

In order to analyze the strong coupling regime (Bq <C 
uj) we apply the canonical transformation (|22|l but now 
with 9 = 2^3. sinwi (and <f> = 0). The new Hamiltonian 
reads 

2B 3 



K = 2B \/l - Q 2 cos P+ — -sincji (55) 



and the equations of motions are given by 



2B 3 



= -2B„^/l - Q 2 sin FH sinujt , (56a) 



P = 2B, 



2B 3 

— cos [ P H sin ujt I . (56b) 



If BqT is small these equations will be in the form 

x = ef(x,u>t), (57) 

with eT< 1, and the averaging theorem can be applied 
|24l |. The mean solution, x m , satisfying \x — x m \ < ce, 
being c independent of e, for t £ [0, =], will be given by 



Xm — &f {Xrn)i 

f°(x) = ±ff(x,ujt)dt, 



(58a) 
(58b) 



where T = ^ is the period defined by f(x, ujt + ujT) = 
f(x,ojt). 

By applying this scheme to the equations of motion 
(|TjFj|) we get 



= -2uj ^1 - Q m 2 sin P m , (59a) 



P rn — 2 ujq 



sfl - Qr, 



: COS P r , 



(59b) 
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being 



2B 3 



(60) 



Jo(x) is the Bessel function of first kind. Notice that we 
may have dynamical localization (l oq = 0) for values of 
S3 and u! that produce zeros in Jo (25L l26| . 

Equations above were obtained by expanding the 
(co)sine in terms of Bessel functions sums. The system 
(|59|l can be solved analytically by noticing that 



fC m = 2ui \J\ - Q m 2 COsP n 



(61) 



is a constant of motion. Note that this result can be 
obtained by applying the formula (|58b|) to Hamiltonian 
(|55(l . Then, since by the averaging theorem Q(t) ss Q m , 
P(t) w P m and JC(t) f=a /C m for t < we can return to 
the original set of variables to obtain the analytical solu- 
tion. However, for our purpose, its sufficient to establish 
the connection between the original and the new Hamil- 
tonian at the stroboscopic times. For all t satisfying the 
condition of validity of the averaging theorem, namely 
B Q T = < i j we have K m = H+ 2B 3 qcosujt, which 
allows us to conclude that 



T~tk = K-m — 2B 3 q k . 



(62) 



This equation allows us to obtain the analytical expres- 
sion for the stroboscopic map in the strong coupling 
regime: 



2^oV^ 



q k * cosp fc , 



(63) 



which could had been obtained from (|61|l by noticing that 
(P m (t k ),P m {t k )) = (P k ,Q k ) = (p k ,q k )- 

Direct comparison with 1)32(1 and 1)52(1 indicates that 
this regime is equivalent to that of high frequency, since 
here we have 7 = 2B 3 and (/) = 0. Then, once more 
the map will be completely symmetric in q and p, with 
central periodic orbits, as those of the Fig. n-(e,f). The 
strong coupling regime has been studied in a more general 
case (quasi-periodic fields) in 

It is also worth to observe that the averaging theorem 
must be valid for all times and not only until t max = -g^ . 
If not, the map would mix two types of contours, one 
determined by (/) = and other by a general (/) 7^ 0. 
In this case, tori crossing would be inevitable. 

We emphasize this point by means of the numerical 
calculation shown in Fig. The plot shows the time- 
independent average (/), calculated numerically from the 
whole evolution (valid for all time), as a function of 

JT we 



1 

BnT 



j^-g-- For small values of t r _ 



with the averaging theorem. Since (/) is the value valid 
for all times (and not only until t max ) we may conclude 
that the averaging theorem, proved valid for large times, 
must be valid for all times. 
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FIG. 7: The average weighted up by the potential time-rate, 
(/}, as a function of = -gKp = the maximal 

normalized instant for the which is valid the average theorem 
solution. In this simulation we fixed B3 = 1.0, u = 10 and 
made vary Bq. According to the averaging theorem (/) = 0. 



3. Resonant Weak Coupling Regime: RWA 

In this section wc apply the averaging theorem to treat 
the NR system in a regime defined by both a resonance 
condition and one small parameter. These two assump- 
tions are essential to the application of the RWA. Ac- 
tually our aim is showing that the averaging theorem 
ensures the validity of the RWA. 

Consider the Scrodinger equation corresponding to 
B^vfl. After a time independent rotation of tt/2 around 
the y-axis we obtain 



~dt 



2(B a z ~ B 3 cosu>ta x )ip 



(64) 



We now perform a unitary transformation to the rotating 
frame, 



have (/) ^ 0, but at such regime the averaging theorem 
is no longer applicable. In fact, according to the numeri- 
cal results (/) already differs from zero for t max /T w 1.6, 
which corresponds to BqT sa 0.6, whereas the theorem is 
valid only for B () T <C 1. Consider now the region of large 
values of t max /T. We see that (/) is null, in agreement cos uitR a x R 



with 

obtaining 

dip' 
l ~dT ~ 

Now, 



ip(t) = R(t)iP'(t), 
R(t) = c^p{-tu>ta z /2), 

(T z — B3 cosiut \R~ l o- x R 



(65a) 
(65b) 

i/>'.(66) 

a x cos 2uit — <T y sin 2uit ) . (67) 
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The second term at the r.h.s. of (|67|l produces the high 
frequency oscillations. According to the heuristics |27j . 
they may be neglected for the purposes of ascertaining 
the "mean" behaviour in time. We now show the condi- 
tions under which this is true. Assume that 



2B n 



(resonance) 



and 



(68a) 



(68b) 



Bs/u) -c 1 (weak coupling). 

Putting (E3 into JSS} we obtain 

dip' S3 / \ . 

i— — = — — ya x + a x cos 2ujt — <j y sin 2ujt J ip . (69) 

Under assumption i|68b[) . it follows, from the averaging 
theorem (Theorem 11.2, pg. 154): 



W(t)-im\=0(B a /u;) 



(70) 



in the time-scale Oitu/B^), i.e., for all t G 0, rgfr^ 
where c is a constant, || • || is the Euclidean norm in spinor 
space and V'oW satisfies the "averaged equation" 



#0 B 3 



From |63]l. (fTHf) and JHJ wc see that 

Mi) ~ i>RWA(t)\\ = O (B 3 /w) , 



for i G 0, 



S37: 



where 



= R{t) e 



lB 3 tcr x /2 



^(0). 



(71) 



(72) 



(73) 



It is possible to estimate the effect of being slightly off 
resonance by (24| (Theorem 11.2, pg. 154). By and 
(|67|l . defining the Rabi frequency 



tt R = [{2B -uj) 2 + S3 2 ] 



1/2 



(74) 



then ip' given by l|6"5|) is close (in the sense of J72J) to the 
solution 



-iQ R t it/2 



where 



V>'(o), 



_ (2S -g) 



(75) 



(76) 



with er 2 = 1. It is to be remarked that an exact method 
suggested by the RWA was devised in [2^. This power- 
ful method (called in j2^| rotating wave transformation) 
goes well beyond the simple idea of the RWA, and has 
recently been applied to the quantum case (|16|) in [l4j . 
We have been concerned, in this section, with the stan- 
dard problem of justifying the replacement of B^, given 
by H20|) . by Bjvit, given by J3UJ), for which, surprisingly, 
no precise estimates on the parameters seen to have been 
given. 



VI. UNITARY QUANTUM NOT OPERATION 

In this section we focus on the potential capability of 
the two-level systems in working as quantum NOT gates. 
We are interested mainly in determining the conditions 
for the realization of the NOT operation, pointing out the 
alternative regimes for its experimental implementation. 

General principles of quantum computation req uire 
that NOT operation be unitary [2^ . References (2iJ |3(j 
propose different constructions of the unitary NOT op- 
eration, e.g., working in an extended Hilbert space. Al- 
though these formalisms are also appicable to our case, 
we are, primarily, concerned with the dynamics We thus 
choose to keep the original setting, in particular, the same 
Hilbert space. In this case, unitarity of the NOT opera- 
tion may restrict us to special cases of initial superposi- 
tions. It indeed occurs, as we will see. 

Our analysis starts by the analytical study of the rotat- 
ing field treated in section llVl for the which all analytical 
solutions are available (appendix [SJ . Since a resonant 
NR field in the regime of weak coupling can be mapped 
on a system governed by a rotating field, as showed in 
the last section, this case will be automatically included 
in the analysis. 

The unitary quantum NOT operation is characterized 
by the transformation S(0) — > S(t not ) in the classical 
phase space. This transition, achieved by means of a 
unitary Hamiltonian dynamics, leads the initial vector 
S(0) to its antipodc S(t not ) in the unit sphere. t no t, de- 
fined by S(t no t) ■ S(0) = S no t • So = — 1, is the instant at 
which the NOT operation occurs. 

Thus, imposing S not • So = — 1 in the analytical result 
given in IA6I we find several regimes of operation for the 
NOT-gate. Bellow we list all of them. 

CASE 1 

Satisfied the resonance condition 



(77a) 



which corresponds to w = -j; the NOT operation will 
occur periodically at 



(2n+l)- 



(n G 



(77b) 



for the following (equivalent) sets of initial conditions 



(PO)?o) 

So 

IV>o) 



(V,0), 

(cosp ,sinpo,0), 
^(l + > + «1- 



(77c) 
(77d) 

(77e) 



Note that the only participants of the unitary NOT op- 
eration are initial conditions belonging to the plane xy, 
specifically on the equatorial line of the unit sphere. This 
situation corresponds to a kind of universal phase NOT- 
gate, since all relative phase po will be transformed into 
Po + tt at t not . We will refer to these solutions as "6<p" , 
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denoting the branch in which all if is valid (equatorial 
line). From now on we will use the notation "fepo" to de- 
note the branch of solutions corresponding to the initial 
conditions related to po- 



CASE 2 

Under the resonance conditions 



B, 



3 — 7-7 



Bo = (2m +1)- (meN), 
the NOT operation occurs periodically at 
tt\ = (2n + l)- („GN) 

UJ 

for the following initial conditions 



(pa, go) = (V,Tl), 
S = (0,0, ±1) 

llto> = l±>, 



(poles), 



and 



(po, 9 o) = (iTr.V) (4 = 0,1,2), 



S = ((-l/v^^^^go), 

l^o) = 



(78a) 
(78b) 

(78c) 



(78d) 
(78e) 
(78f) 



(78g) 
(78h) 



^l+) + (-l)S/^|->. (781) 



The set of initial conditions referent to the poles has 
already been reported in Q • Notice that the points corre- 
sponding to the poles in the unit sphere are transformed 
in separatrices (q = ±1) in the phase space qp. It is just 
a manifestation in pq of the ambiguity related to <p in the 
poles of the unit sphere. 

The last set corresponds to the case of a universal real 
NOT-gate, since the initial quantum state is composed 
only by real coefficients. Note that there exist two dif- 
ferent branches, depending on the value that I assume. 
The branch 60, for I = 0, 2n, contains all initial condition 
satisfying Si > 0. On the other hand, for the branch bir 
(I = 1), one has Si < 0. These branches meet each other 
at Si = (<7o = ±1, poles), onto which pass a separatrix. 

CASE 3 

Provided 



B, 
B 



UJ 

2' 
w, 



the NOT operation occurs periodically for 



An) 
L not 

(Po, Qo) 
So 



(2n+l)— (net 
(V,Tl), 

(0,0, ±1) (poles), 
l±>, 



(79a) 
(79b) 



(79c) 

(79d) 
(79e) 
(79f) 



and for 

An) . 

°not ~ 

(po,qo) = 

So = 

|V>o) = 

in which 



(4n + e t 
[2l + 3]j,V 



(n G 



(79g) 
(79h) 



1 - So 1 



(V 1 - 9o 2 cosp , V 1 ~ 9o 2 sinp ,-9o) , (79i) 

-), (79j) 



(79k) 



1 + qo 



1, for / even, 
3, for I odd. 



Note that in this case, t no t depends on the choice for the 
Po- Once again, there are several branches contributing 
to the NOT operation, namely, 6f , b^-, b^- and b^-, 
each one of them corresponding to a given I. These cases 
correspond to quantum states with complex phases in the 
form e' lp " = (±1 ± «)/v2, for all combinations of signals, 
for any qo- Therefore, also in these cases, there exist a 
type of universality for the NOT operation. 



CASE 4 

With 



2' 

UJ 

4m 



B 3 = 
Bo = - 



the NOT operation occurs periodically for 



An) 

L not 

So 

IV'o) 



n 2tt 

m(2n + 1) — 

UJ 



(n G N), 



(v/l-go 2 ,0,- go ), 



1 - qo 



\+)+* 



1 + 90, 



(80a) 
(80b) 

(80c) 
(80d) 
(80e) 
(80f) 



/AjT 



This case corresponds to the branches 6-| and b- 2 . 

For completeness, we point out that all informations 
got above concerning the initial conditions emerged from 
the following simple expressions derived from l|A6|l at a 
specific t not : 



(1) S not -S = l-2q, 



(81a) 
(81b) 

qo 2 + (qa 2 - 1) sin (2p ) , (81c) 

(81d) 



(2) S not ■ S = -qo 2 + (qo 2 - 1) cos (2p 

(3) Snot ■ So 

(4) Snot ■ So 



-qo +(l-qo ) cos (2p ) 



in which we have to impose that S no t ■ So = —1. Figure 
IH1 presents an overview of the global symmetry presented 
by the set initial conditions allowed to implement the 
unitary NOT operation. 

It is interesting to note that the initial conditions found 
in the case 4 are just localized over the separatrix of mo- 
tion in phase space pq. However this is not relevant to 
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FIG. 8: Set of initial conditions able that allows for the imple- 
mentation of the NOT operation in the R system (solid lines 
and the points corresponding to the poles of the sphere). The 
solutions are shown in both the phase space pq (on the left) 
and the tridimensional phase space xyz. By simplicity, only 
the north hemisphere is shown, since the south presents the 
same contents. 

the quantum dynamics, since the separatrix exists only 
in this particular phase space. It occurs in virtue of 
the topology change from a tridimensional phase space 
(sphere) to a bidimcnsional one (plane qp). A deeper 
analysis of the classical picture associated to the quan- 
tum dynamics may be found in appendix FBI 



A. NOT operation in the NR system 

The precedent analysis shows that practically all NOT 
regimes occurs for the resonance B 3 = j , which produces 
stroboscopic maps with centered periodic orbits. The 
connection between R and NR maps verified in section 
IVl indicates that we have to search regimes such that 
£?3 = 5 in the NR system, since the parameter 7 plays 
the role of the R frequency lo for the map symmetry 
(see equations and I33[) . According to (|51|l and i|52|) 
this regime will be achieved only asymptotically, provided 
Lo 2 (B a 2 + i?3 2 ) -1 —> 0. We have verified numerically 
for the case 2 (Fig. [§J| that the NOT operation indeed 
occurs asymptotically, i.e., S rl0 4So — ► —1 as — ^ — * 0, with 
B3 = B . 

The most interesting and non-trivial regime of NR 
NOT operation was found at a resonance equivalent to 
the case 1. Obeying the analogy between 7 and lo we 
obtain the following resonance condition for the NR sys- 
tem: 

7 2 = Bo 2 + (B 3 - |) 2 , (82) 

being 7 a parameter determined by the choices of the 
NR frequencies lo, B3 and Bq, i.e., 7 = 7(0;, Bq, B3). 
We determined numerically the parameters satisfying 
(|82|l by fixing lo and B3 and varying Bq. The result 
is shown in Fig. ^] Once found the solutions for equa- 
tion (|82J) , we verified the realization of the NOT opera- 
tion for the NR system for a universality class given by 
(Po,Qo) = (R + |] Vj W1 th I G N, as can be seen in 
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FIG. 9: NR NOT operation with B = B 3 = 0.2 and (a) 
lo — 3.0 or (b) lo = 10.0. 7 tends to 2B3 asymptotically, thus 
satisfying the resonance conditions of the case 2. However, 
here the universality class is (po,9o) = (n?>V). The vertical 
dashed line stands for the analytical NOT instant given by 
tnot = 7.854 (see appendix IBl. 
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FIG. 10: Numerical determination of the Bo value that satis- 
fies 1821 (tick solid line) for lo = 1.0 and B3 = 1.5. For each 
value of Bo we determine 7 (bullets) by means of the fitting 
procedure described in section W\ 

Fig. ITTKa). We also plot the R NOT operation with 
analog parameters. 

Although not shown by the figure, the NR NOT oper- 
ation is also periodic, with t no t = (2n+l)57r. Notice that 
the first occurrence of the NOT operation happens in the 
NR system five times latter than in the R system. This 
is an adverse point for the NR system, since decoherence 
effects may be not negligible in this time scale. 

At last we point out that the prediction of NOT op- 
erations in the NR system was based on the comparison 
of stroboscopic maps with those of the R system. These 
classical mathematical tools, based on the concept of tra- 
jectories flowing in the phase space, were crucial in the 
determination of the connection between the R frequency 
lo and its NR counterpart 7, parameters that define the 
symmetry of the maps. It stresses the usefulness of the 
classical gyromagnctic picture, being an example of clas- 
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FIG. 11: (a) NR NOT operation for u) = 1.0, 7 = 1.486, 
B3 = 1.5 and Bo = 1.279 occurring at t not = 5tt for several 
initial conditions with po = §■ (b) R NOT operation for 
u = 1.486, B 3 = 1.5 and B = 1.279, with t not = (2n + 1)£ 
for initial conditions given by (po,9o) — (V, 0) (case 1). 

sical analysis allowing quantum predictions. 

VII. SUMMARY AND FINAL REMARKS 

We have presented some relevant results concerning 
the underlying classical aspects referent to the two-level 
systems under time-dependent external fields. The main 
results were obtained in virtue of the application of the 
classical gyromagnet picture, which has been shown to 
be a powerfull approach to investigate quantum two-level 
dynamics. 

We have shown that the underlying classical dynamics 
is intcgrable for arbitrary external fields as a consequence 
of the unitary quantum dynamics. The numerical stro- 



boscope maps, which are in perfect agreement with our 
analytical proof of integrability, were crucial mainly to 
the identification of the similarities between the R and 
NR dynamics. On the other hand, the average weighed 
up by the potential and the application of the averaging 
method allowed us to identify several interesting regimes 
with explicit analytical solutions. 

It also should be remarked that we have arrived at 
other two results of independent interest: (i) it has 
been established beyond doubt (although not by rigorous 
proofs) that both the RWA and the averaging theorem 
are valid for such systems well outside their expected re- 
gions of validity; (ii) we have shown that the RWA is a 
formal consequence of the averaging theorem. 

Concerning the realization of the NOT-gate in two- 
level systems, we have found several regimes for its ex- 
perimental implementation for two types of external pe- 
riodic fields. For a rotating field, the quantum dynamics 
has been solved and universality classes for the initial 
quantum state have been identified analytically, thus in- 
dicating all control parameters for the NOT-gate work- 
ing. On the other hand, for the nonrotating field, analyt- 
ical equations of motion were not found and predictions 
about the NOT operations were possible thanks to the 
classical analysis based on stroboscopic maps. 
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APPENDIX A: ANALYTICAL SOLUTIONS 

As we mentioned, the solutions for the rotating field 
treated in section IIVI were found by solving the quan- 
tum problem. Applying usual unitary transformations, 
namely exp (— luita^), we led the system to the rotating 
frame, in which the system is no longer time-dependent. 
Then the solution was promptly obtained. The last step 
consisted in applying the parametrization for classical 
canonical variables, S = {S\, S2, S3) and (p, q). The fi- 
nal solutions, including that one for S(t) • S(0), are given 
bellow as functions of the initial condition (po,qo), with 
B = 2 v / iV r +TF and Q = B 3 -%. 
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Si(t) 



S 2 (t) 



s 3 (t) 



4 go Bp 
B 2 



9 B 

2 il cos(u) t) s'm(B t/2) + — sin(iu t) sm{B t) 



B 2 



2 B 2 cos(p + wt) + 



— - Q j cos(p + wt-Bt) + + C0S (P° + wt + Bt)+4 B 2 cos(p - tut) sin(B i/2) 2 



Al) 



4 go 5q 
B 2 

B 



9 B 

217 cos(u)t)sin(Bt/2) - — sin(tyt) 8111(5 1) 



B 



B 2 



2B 2 sm{p Q +wt) + 



tt) sin(p + tut - Bt) + I y + ft ) sin(p + tut + Bt) - 4_B 2 sin(p - wt)sin(Bt/2)' 



(A2) 



^ [ft 2 + cos(Bt)B 2 ] + 4g o\/l n cosOo) [l- cos^t)] +| sin(po) sin(St)}. (A3) 



Pit) 



S 3 (t). 



arctan 



S2(t) 

[Si(t) 



0, if5i(t)>0; 
7T, ifS"i(t)<0. 



(A4) 
(A5) 



S(t) • S(0) 



\B 2 cos(tut) + 2 + ft 2 ) cos(wt) cos(Bt) -2Btt sin(tut) sin(Bi) 



'go 



B 2 

cos(B t/2) sin(f) + ft cos(^) sin(B t/2) 



n 2 



B 2 



+ 



8B q y/l - q 2 cos(p - 2 cos(i^) sin(B t/2) 2 + f sin(^) sin^t) 



B 2 



4B, 



o 2 sin(B t/2) 2 |cos(2p — tu t) — q 2 1 + cos(2p — tu t) j 



B 2 



r 



(A6) 



APPENDIX B: GEOMETRIC PICTURE 

The classical version of to the two-level quantum sys- 
tem may be understood as a classical gyromagnet, with S 
precessing around the external field B. In this appendix 
this analogy is emphasized by means of a geometrical 
picture of the classical dynamics, which illustrates ques- 
tions concerning the separatrix, the periodic orbits and 
the NOT operation. 

We start our analysis by focusing on the R system in 
its rotated version, namely, 



JC = — B • S = —B cos*, 

B = (B x ,B y ,B z ) = -2 (5 , 0, ft), 



(Bla) 
(Bib) 

S = Wl - Q 2 cos P, \/i - Q 2 sin P, -Q), (Blc) 

with B = ||B|| and being * the angle between the exter- 
nal field and the unit vector S. The dynamics is governed 
by 

v,f = SxB, (B2) 



which defines the velocity of the unit vector S. We may 
calculate the acceleration a = v noting that the external 
field is timc-indcpcndcnt. This yields 



a = v x B. 



(B3) 



Consider two orthogonal versors, say en (parallel to the 
B) and ej_ (perpendicular to B), such that we may de- 
compose: S = cos* e | + sin*e^ and B = -Ben. Since 
JC and B are constants the only option for S is perform- 
ing a precession around the static magnetic field. This 
consideration makes the versor time-dependent. 

Consequently, the velocity and acceleration can be 
written respectively as 



v = B sin * (e^ x e 
a = B 2 sin* (-ej_) . 



\\) ' 



(B4a) 
(B4b) 



Being JC and B constants, we conclude by (|Bla|) that *, 
v = ||v|| and a = ||a|| are all constants too. Moreover, 
(|B4b|l indicates that a is centripetal, stressing the fact 
that S develops a precession motion around B. 
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By (|Bla|) wc note that * G [0,tt] and K G [~B,B]. 
This allows us to re-write equations l|B4[) in terms of the 
energy as 



\/B 2 - JC 2 (e ± x e 



) ■ 



a = B\J B 2 - K 2 (-ej_). 



(B5a) 
(B5b) 



Finally, we may define the angular velocity w = 
(sin "J ei) x v which defines the precession direction. Us- 
ing equations given above we get 



Bsin 2 * 



\B — JC 



II- 



(B6) 



which shows that the angular velocity is always anti- 
parallel to B. This helps us to understand the direction 
of the trajectories flux. Moreover, since w = ||w|| is con- 
stant we may define a rotation angle as = <& + 



1. Periodic Orbits (Eigenstates) 

A first interesting consequence of the simple scheme 
developed up to here is that the initial condition corre- 
sponding to the periodic orbits (eigenstates) of the sys- 
tem may be promptly obtained for K? = B 2 , thus yield- 
ing v = a = 0. Therefore, these periodic orbits corre- 
spond to "static" quantum states (eigenstates). In this 
situation, S and B are parallel = 0,7r), and there is 
no external torque. 



2. Separatrix 

The classical gyromagnet presents in the space pq a 
separatrix of motion associated to the specific energy 
±2fi. Therefore, by HI with Q = B 3 - f and = 0, the 
explicit contour for the separatrix at stroboscopic instant 
is given by 



Bq^/1 - q k 2 coap k = (±1 + q k ) ft. 



(B7) 



This equation offers all forms for the separatrix, as is 
shown in Fig. 1121 



q 




IC.dL 



q 




q 



2p 



2p 



FIG. 12: Separatrices of motion plotted with the same pa- 
rameters of the respective graphics shown in Fig. 

We then have JC = ±B Z , in which the signals indicate 
the direction of the vector S compared to B fixed. Con- 
sider now that O is the projection angle of the magnetic 



field along the direction z, i.e., B z = BcosQ. Hence, ac- 
cording to (|Bla() we obtain ^ , = 6or5' = 7r — O, which 
indicates that the vector S will be necessarily parallel (or 
anti-parallel) to the versor z at some instant. Therefore, 
inevitably S will pass through the poles during its pre- 
cession motion. However, concerning the dynamics of the 
unit vector S on the sphere, there is absolutely nothing 
special in this situation. For instance, the period of the 
precession in this case is 

2lT 2-7T 7T 

r = ^ = s7 = W (B8) 

which is remarkably smaller than the period (infinity) ex- 
pected for a separatrix. In this sense, there does not exist 
a separatrix in the classical unit sphere. It appears only 
in the pq (or PQ) phase space reflecting the peculiarities 
of its topology, which transforms "localized" states on 
the sphere (poles) into "delocalized" ones on the plane 
(lines (p,g) = (V,±l)). 

In addition, we note that the hyperbolic fixed points 
(^, ±1) by the which the separatrix passes through, have 
no special physical interpretation on the sphere. Figure 
II 31 illustrates geometrically our argumentation. 




FIG. 13: Unit sphere projected in the plane xz. For the 
separatrix energy K, — ±B Z the vector S precesses around B 
with * = O. In this situation the unit vector crosses the poles 
of the sphere, which corresponds to the unstable points of the 
space PQ. The angular velocity w and the flux direction are 
also indicated. 



3. NOT operation 

Let us work with matrix notation, using 5^ = 
(S x S y S z ) and Sj' = (Si S2 S3) for the line matrices 
corresponding to the unit vectors in phase spaces PQ 
(index R) and pq (index r) respectively. The superscript 
T denotes the transposition operation. According to the 
canonical transformation (|22|l we may write 

S R (t) = G(t)S r (t), (B9a) 
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(cos Lot sin tot \ 
-sin Lot cos Lot , (B9b) 
ij 

being G(t) a rotation matrix around the axes z (or "3"). 
Consequently, G -1 (t) = G T (t) = G(-t) and S r (t) = 
G(—t)S]i(t). Now, it is easy to see that 

Sl(t)S R (0) = S?(t)G(-t)S r (0), (BlOa) 
5 r T (t)5 r (0) = Sl(t)G(t)S R (0). (BlOb) 

This result shows that in general the NOT operation does 
not occur simultaneously in both frame. 

The NOT operation instant T not in the frame R may 
me determined from the geometrical picture developed in 
this appendix. It reads r not = S = 5 j which is a half of 
the precession period with \V = ^. This value for <3> is a 
necessary condition for the NOT operation in the frame 
R. 

All regimes found for the NOT operation in section IVT1 
may be formulated in terms of a single relation. To get 
it we re- write Ea. l|Bla|) at t = as 

£cos* = B x Si(0) +5^(0), (Bll) 

in which S*i(0) = S x (0) and 5 3 (0) = 5 Z (0), since 
(-PcbQo) = (PotQo)- br au regimes identified we have ei- 
ther 53(0) = (q = 0) or B z = (resonance B 3 =u>/2). 
Then, using the projection angle O defined in section lB~2l 
we get 

cos* = Si (0) sin 9, (B12) 

with 5i(0) = y/T — <7o 2 cospo- Relation (|B 12|l puts to- 
gether informations about the precession angle * (related 
to the energy), the magnitudes of the magnetic field pro- 
jections (related to resonances) and the initial conditions 
(related to the initial quantum state) . All regimes of the 
NOT operation obey this simple formula. 

By (|B12|) we see that the NOT operation will occur in 
both frames either if (po,qo) = (f)V) or if (po,<7o) = 
(V, ±1), since B and O are always different of zero. 
This conditions are verified in cases 2, 3 and 4. Fur- 
thermore, in theses cases B3 = lo/2, so B = 2Bq and 



Tnot = 2S - = tnot- Therefore, the NOT operation will 
occur simultaneously in both frames. 



4. Mean NR. NOT operation 

Let us consider the NR system, for the which the mag- 
netic field is given by B(t) = — 2(B , 0, B 3 cosLot). In the 
high frequency limit we may apply the averaging theorem 
to get 

U = B • S (B13a) 



T 

B = l jB{t)dt={-2B ,0,0). (B13b) 


We see that in average, the field behaves like a time- 
independent external perturbation, allowing us to use the 
gyromagnct picture developed above. Thus, we may de- 
scribe the precession motion by means of a constant av- 
erage angle * such that H = —Bcos^>, with B = ||B||. 
This equation may be written at the initial instant with 
(|BT3|) as 



COS* = —y/l - q 2 COSpo- (R14) 

Following the steps of the precedent sections one may 
determine the angular velocity. The result is simply 
w = Bsin*. But we know that the NOT opera- 
tion occurs only with \V = ^. In this situation, the 
NOT instant is directly calculated from the expression 
for w as the time spent for performing a half rotation: 
tnot = ^ = Moreover, in this regime equation (|B14|) 

reads yl — qo 2 cospo = 0, which is satisfied in general 
for {po,qo) = (R + 5] TTjV). This result explain the uni- 
versality class found in section IVll for the NR system in 
the high frequency limit. Furthermore, using the numer- 
ical parameters of the Fig. [5] we get t not = 7.854, in 
perfect agreement with the numerical results. 
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